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ABSTRACT 



Aims. We present an analysis of simultaneous X-Ray and UV observations of comet C/2007 N3 (Lulin) taken on three days between January 
2009 and March 2009 using the Swift observatory. 

Metiiods. For our X-ray observations, we used basic transforms to account for the movement of the comet to allow the combination of all 
available data to produce an exposure-corrected image. We fit a simple model to the extracted spectrum and measured an X-ray flux of 
4.3±1.3 X 10"" ergscm"^ s"' in the 0.3 to 1 keV band. In the UV, we acquired large-aperture photometry and used a coma model to derive 
water production rates given assumptions regarding the distribution of water and its dissociation into OH molecules about the comet's nucleus. 
Results. We compare and discuss the X-ray and UV morphology of the comet. We show that the peak of the cometaiy X-ray emission is offset 
sunward of the UV peak emission, assumed to be the nucleus, by approximately 35,000 km. The offset observed, the shape of X-ray emission 
and the decrease of the X-ray emission comet-side of the peak, suggested that the comet was indeed coUisionally thick to charge exchange, as 
expected from our measurements of the comet's water production rate (6-8 xlO^** mol. s"'). The X-ray spectrum is consistent with solar wind 
charge exchange emission, and the comet most likely interacted with a solar wind depleted of very highly ionised oxygen. We show that the 
measured X-ray lightcurve can be very well explained by variations in the comet's gas production rates, the observing geometry and variations 
in the solar wind flux. 
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1. Introduction 



Comets emit X-rays via the process of solar wind charge ex- 
change (SWCX). Gas in the coma of the comet donates one or 
more electrons into an excited energy level of a highly-charged 
ion of the solar wind. In the subse quent relaxa tion of the ion 
a UV or X-ray photon is emitted (ILisse et al 
1 997t iKrasnopolskvll 1 997l) . 

SWCX in comets probe s states of the 
throughout the heliosphere dSchwadron & Cravens 



This sample contains a broad variety of comets, solar wind en- 
vironments and observational conditions and clearly demon- 
strates the diagnostics available from cometary charge ex- 
change emission. Several of those comets showed interesting 
large scale s tructu res in X-ray such as the jets in 2P/Encke 
(ILisse et al.l 12005), the Deep Impact triggered plume in 



— I — — — 1 (ILisse et at.l IZUUJ | ). tne Deep impact tnggerea plume in 
1996t i CravensI gp/jempel 1 jLisse et al1l2007h . and the disintegrating comets 



solar wind 
2000[ 



73P/Schwassmann -Wachmann 3 and C/1999 S4 (LINEAR) 
dWolk et al.1 12009|) . These interpretations relied on ground- 



Kharchenko & Dalgamd 200 ih. To date, over 20 comets have 



based observations to initially identify those structures. Several 



been observed in X-ravs (ILisse et al . | 1996|: iDennerl et al 119971: 
Lisse etalj|2004t iKrasnopolskvl 120061: iBodewits et alj|2007h 



Send offprint requests to: J.A. Carter 



200C; Lisse et al. 2005 



(Lisse et aL 1999: Neusebaueret al. 


Willineale et al. 


20061 


Lisse et al. 



20071) via relationships between the comet gas production rate. 
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the heliocentric distance, and the behaviour of the solar wind. 
This paper however presents the first endeavour to directly con- 
strain the cometary gas production rates and relate this to the 
observed X-ray variability by employing Swift's co-aligned in- 
strument suite, allowing us to make simultaneous X-ray and 
UV measurements from a single spacecraft. 

C/2007 N3 (Lulin) is a dynamically new comet that 
was discovered by Lin Chi-Sheng and Ye Quanzhi at Lulin 
Observatory, Taiwan, in 2007. In this paper we analys e data 



pie are indicators of ICME plasma dRichardson & Cane 2004 



from comet Lulin obtained by the Swift satellite (Gehre ls et al 



|2004 on several days in early 2009 as part of its Targets of 
Opportunity programme for non-gamma ray burst targets. The 
days were chosen to sample emission from the comet at various 
stages during its passage into and out of the inner solar system, 
taking advantage of its high activity and its close proximity to 
Earth. An initial analysis of these ob servations using data fr om 
the Swift UV grism was discussed in lBodewits et al.l (1201 ih . 

Comet Lulin moves in an orbit of low inclination of just 
L6° from the ecliptic, allowing us to link measured solar wind 
data from the Solar Terrestrial Relations Observatory space- 
craft (sTmEOj_spacecraft A, part of the plastic plasma pack- 
age, |BlushetalJ|2005[) to our observations. We expected that 
this comet would sample the low-latitude, highly-charged so- 
lar wind and that the intensity of this emission would reflect 
changes in the solar wind flux, the quantity of neutral species in 
the comet's atmosphere and the distance to the observer, sim- 
ilar to wha t has been seen i n other near-ecliptic c omets (e.g. 



2P/En cke (iLisse et all 119991 120051) : P/Tempel 1 (ILisse et al 



20071)). 



The observations were taken around solar minimum, when 
the solar wind can be simplified to the stratification of a slow, 
low-latitude wind originating about the solar equator and a 
higher-latitude, faster wind. The heavy-ion freeze-in temper- 
ature at the solar corona determines the ion composition of 
the out-flowing solar wind. In the bimodal state outside solar 
maximum, the polar wind ha s a lower io n freeze-in state than 
the wind at lower latitudes ( Geiss et al] |l995). The nominal 
low-latitude solar wind flows radially away from the Su n, yet 
the Su n is rotating, resulting in a Parker spiral outflow ( Parkei 



19581) . However, this flow contains streams of plasma with dif- 



ferent radial velocities. Co-rotating Interaction Regions (CIRs) 
occur when a faster stream piles up against the slower, ambi- 
ent plasma, resulting in an upstream compression region and 
a downstream rarefaction region. CIRs are characterised by an 
increase in ion density (with a thickness of ~ 0. 1 AU) followed 
by a decrease in this dens ity, and a correspondin g sharp in- 
crease in the ion velocity ( IGosUng & Pizzo 1 19991 and refer- 
ences therein). They are found at low to mid heliospheric lat- 
itudes. Coronal Mass Ejections (CMEs, or ICMEs when they 
are found within interplanetary space) however, do not form 
part of the nominal co-rotating solar wind flow, but are large 
clouds of fast moving plasma ejected from the Sun with dis- 
tinctive compositional signatures. Along with differences be- 
tween the fast and slow wind compositional fractionation as 
described above, transient events such as ICMEs show marked 
abundance signatures that can be used to identify plasma of 
this type. The presence of highly-charged iron, elevated oxy- 
gen states along with enhanced a-to-proton ratios for exam- 



Zurbuchen & RichardsonlEooet IZhao et al.ll2007l) 

This paper is organised as follows. In Section|2]we discuss 
the data analysis steps taken for both the X-ray and UVOT data 
sets. In Section[3l we discuss our background and spectral anal- 
ysis of the XRT data along with the X-ray radial extent of the 
comet, the morphology of the X-ray and UV emission, the gas 
production rates of the comet and the X-ray temporal variabil- 
ity. We finish with our conclusions in Section|4] 

2. Data reduction 

2.1. X-ray telescope 



We use data from the X-ray telescope (XRT; iBurrows et al, 
|2005) in photon counting mode. The XRT is sensitive over the 
energy ra nge 0.2 to lOkeV , has an effective area of ~80cm^ 
at 1 keV ("Godet et al."2009'), a spectral resolution of ~130eV 
at Si-Ka (Pagani et al. 2011) and a sensitivity limit of 2 x 
10 ''^ergscm"^ s The primary source of the XRT particle- 
induced background is protons incident on the focal plane, 
swept up as Swift moves in its low-Earth orbit. 

The 14 Swift XRT observations used in this analysis, taken 
on three different days (2009-01-28, 2009-02-16 and 2009-03- 
04, 2009 UT) are listed in Table 1 . These observations were at 
fixed pointings and Swift did not track the comet. Ephemeris 
data for the comet were obtained from the Horizons websiteQ. 
We used the cleaned XRT event lists which have been filtered 
for periods of high background by the standard Swift process- 
ing chain. The event files had also been adapted for an improve- 
ment in the gain function of the XRT (XRT instrument team, 
private communication). No bright X-ray astronomical point 
sources were formally detected in the field-of-view (FOV) for 
any of the pointings by the source detection algorithm detect of 
the HEASARC task ximagiH . 

For every event file we transformed each individually 
recorded X-ray photon in detector space, using the ephemeris 
information to account for the comet's movement across 
the sky. We calculated the change in Right Ascension and 
Declination of each event and calculated the new detector 
coordinates accordingly. Every exposure was transformed so 
that the nominal position of the comet at the exposure start 
was found at the centre of the detector, and we accounted for 
the change in Sun-direction by rotating each frame. We then 
stacked these event files together to improve the signal to noise 
and adapted the exposure keywords appropriately. Vignetted 
exposure maps were produced for each observation. Each ex- 
posure map was stretched to account for the movement of the 
comet and shifted and rotated as per the steps applied to the 
event files. All exposure maps were stacked to provide a com- 
bined exposure map to accompany the stacked event file. 

We created a smoothed image of the stacked events and di- 
vided this by the stacked exposure map to produce an exposure - 
corrected image. As there was no significant emission detected 
above 1 keV, we study only photons with E < 1 keV here. We 



http://ssd.jpl.nasa.gov/horizons.cgi 
http://heasarc.nasa.gov/xanadu/ximage/ximage.html 
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1 arcmin 





Fig. 1. Upper panel: a smoothed 0.3 to 1 keV exposure cor- 
rected image of the stacked XRT event Hst following the trans- 
formations described in Section 12.11 The circle (green) and 
annulus (red) refer to regions used for spectral analysis, see 
Section im Lower panel: zoomed-in UV image of the comet, 
from the first UVWl exposure of observation 31332 (see Table 
1), taken with a pointing Right Ascension and Declination of 
(231.44° and -17.89°), with X-ray contours from the upper im- 
age overlaid. Contours start at 1.0 x lO"-' ct s"' arcmin"^ (~ 6.5 
sigma above the background), and increase in increments of 0.4 
X lO"-' ct s"' arcmin"^. A one arcminute bar is shown on each 
panel. 



restrict the lower limit of our analysis to 0.3 keV d ue to the un- 



certa in calibration of the XRT below this energy (iGodet et al 
12009 ). A 0.3 to IkeV exposure-corrected image is shown in 
Figure [1] (upper panel), which has been rotated so that the Sun 
vector is as for the first observation of the data set. This image 
has been smoothed employing a Gaussian filter with FWHM 
4.6 pixels (one pixel equals 2.36 arcseconds side length). 



2.2. UV-Optical telescope 

Along with the X-ray observa tions, the Swift- 



Ultraviolet/Optical Telescope (UVOT; iRoming et alJ 120051) 



acquired 75 observations of comet Lulin using its optical and 
UV broad-band filters. UVOT provides a 17 x 17 arcminute 
FOV, with a plate scale of 1 arcsecond/pixel and a point 
spread function of 2.5 arcseconds FWHM ("Ma son et al 



2004) . Additionally, on 2009-01-28 the UV grism was used 
to acquire low resolut ion (A/ 6 A - 100) spectra of the comet 
dBodewits et al.l |201 ih . UVOT images of comet Lulin using 
the uvwl (Ac 2600A, FWHM 700A) and v-band (A^ 5468A, 
FWHM 750A) filters are discussed later in Section IT21 

Water vapour and its fragment species are the most abun- 
dant volatiles in cometary comae. Co met Lulin produced no 



more than 2-5% C O2 and <2% CO dBodewits et alj 12011 



lOotsubo et al. 2010l) . Given that these molecules have com 
parable charge e xchange cross sections at solar wind ion ve- 
locities as H2O dMawhorter et al.ll2007h . the number of water 
and hydroxyl molecules is a good proxy for the comet's con- 
tribution to the variability of the X-ray lightcurve. We used 
Swift's uvwl and v-band filters to measure the number of OH 
molecules in the FOV, and to estimate Lulin's water production 
rate during our observations. 

We measured the comet's flux in large apertures of 400 
arcseconds radius to get excellent signal to noise, to sample 
a large fraction of the coma, and to allow for comparison with 
the X-ray photometry (Section I3.2l l. We obtained coma pro- 
files by azimuthally averaging the surface brightness. This was 
achieved by converting each UVOT image into polar coordi- 
nates (r, (p) and finding the mean surface brightness at a given 
radial distance (r). By taking the azimuthal-averaged surface 
brightness for every pixel at a given radial distance (r) from the 
optocenter (area of brightest emission, that is not necessarily 
coincident with the position of the nucleus), we constructed a 
2-dimensional azimuthal average image of the comet that was 
used to search for faint structures. 

Although this method assumes the large scale coma to be 
symmetric, it enables the effective filtering out of stars and 
other background objects. The coma fills the entire UVOT 
FOV, but by comparing the signal at the edge of the detector 
with Swift calibration data we estimated a background signal 
of 0.005 counts s"' pixeL' (uvwl) and 0.02 counts s"' pixeL' 
(v-band), corresponding to 1% (uvwl) and 13 % of the total 
signal from the comet. At the comet's optocenter, the count 
rates were high enough to result in some coincidence loss (here 
approximately 7% for uvwl, and <20% for v-band ) of the flux 



in a 5 arcseconds aperture around the optocenter (iPoole et al 



20081). We decided not to correct our photometry for this co- 



incidence loss, as it only affected a small fraction of the total 
flux in the 400 arcsecond aperture. Based on these considera- 
tions we estimate stochastical errors to be 5% (uvwl) and 15% 
(v-band). Absolute er rors in the photometry are better than 3% 
ffireeveld et al.ll201 lb . 

The measured fluxes in the filters relevant to this study 
(uvwl and v-band) are summarised in Table 1. Note that two 
UVOT exposures with the same observation ID (31332) were 
obtained in the same Swift orbit, while in observations 31336 



Table 1. Swift XRT and UVOT observations of comet Lulin performed in early 2009, listed by the observation identifier (OBS JD). Start and stop times are given for the XRT exposures in the 
format MM-DD. We also note the XRT livetimes (T), the Sun-comet distance (rn), the Earth-comet distance (A), the Sun-Earth-comet (S-E-C) angle and the rate of change in Right Ascension 
and Declination (ARA, ADec). The horizontal lines separate the three observing periods. We list the UVOT measured fluxes for the uvwl and v-band filters in 400 arcsecond apertures (see 
Sec. 2.2). Flux (OH) is the flux in the uvwl filter after removal for the continuum contribution (see Section 3.3). N(OH) is the number of OH molecules in a 400 arcsecond aperture and Q(OH) is 
the production rate of OH molecules (see Section 3.3.1). No v-band observations were obtained on 2009-01-28, and we list upper limits for the OH production rates based on the total measured 
uvwl flux. Only stochastic errors are given. The systematic uncertainty in the UVOT measured fluxes is less than 5%, and approximately 25% for Q(OH). We list the total number of X-ray 
photons for each pointing (Nx) within the spectral extraction region and between 0.3 and 1 keV, as described in Section 3.1. We also list the measured XRT X-ray fluxes (between 0.3 and 1 keV 
in cgs units (ergs cm"^ s"' ), see Section 3.3.3) in the last column, or upper limits where appropriate (assuming a negligible background). 



OBSJD 


Start 


Stop 


T 




A 


S-E-C 


ARA. 


ADec. 


Flux uvw 1 


Flux v-band 


Flux (OH) 


N(OH) 


Q(OH) 




XRT flux 








(s) 


(AU) 


(AU) 


(deg) 


(sec/s) 


(sec/s) 


(Wm-^sr-') 


(Wm-^sr-i) 


(Wm-^sr-i) 


(molec.) 

(Xl033) 


(molec. s"') 
(xlO^S) 




cgs 
(xlO-") 


31332 


01-28.00 


01-28.08 


2335 


1.245 


1.083 


74.43 


-0.025 


0.006 


0.30 ± 0.01 






<11.0 


<9.65 


273 


10.4 + 3 


31332" 


















0.29 ± 0.01 






<11.0 


<9.63 






31336 


01-28.13 


01-28.14 


1136 


1.246 


1.079 


74.64 


-0.026 


0.006 












138 


<25 


31337 


01-28.19 


01-28.27 


2276 


1.246 


1.077 


74.75 


-0.026 


0.006 












299 


11.0±3 


31347 


02-16.53 


02-16.55 


943 


1.338 


0.508 


125.04 


-0.132 


0.051 


0.59 ± 0.03 


0.81 + 0.12 


0.47 ± 0.06 


3.16 ±0.4 


5.5 ± 0.7 


58 


<6.1 


31348 


02-16.61 


02-16.62 


932 


1.339 


0.506 


125.36 


-0.133 


0.051 


0.60 ± 0.03 


0.86 ±0.13 


0.47 ± 0.06 


3.11 ± 0.4 


5.5 ±0.6 


44 


<6.2 


31349 


02-16.67 


02-16.69 


948 


1.339 


0.505 


125.66 


-0.134 


0.052 


0.60 ± 0.03 


0.85 + 0.13 


0.47 ± 0.06 


3.11 ±0.4 


5.5 ± 0.6 


52 


<6.5 


31350 


02-16.74 


02-16.75 


935 


1.340 


0.503 


125.95 


-0.134 


0.052 


0.60 ± 0.03 


0.85 ±0.13 


0.47 ± 0.06 


3.11 ±0.4 


5.5 ± 0.7 


53 


<5.8 


31351 


02-16.81 


02-16.82 


937 


1.340 


0.502 


126.25 


-0.135 


0.052 


0.60 + 0.03 


0.84 ± 0.13 


0.47 ± 0.06 


3.10 ± 0.4 


5.5 ±0.7 


69 


<9.8 


31352 


02-16.88 


02-16.89 


940 


1.341 


0.500 


126.56 


-0.136 


0.053 


0.59 ± 0.03 


0.86 ±0.13 


0.46 ± 0.06 


3.04 ± 0.4 


5.5 ± 0.6 


82 


<10.2 


31353 


02-16.94 


02-16.95 


932 


1.341 


0.499 


126.86 


-0.137 


0.053 


0.60 ± 0.03 


0.86 + 0.13 


0.47 ± 0.06 


3.05 + 0.4 


5.5 + 0.6 


87 


7.1 ±3 


31364 


03-04.31 


03-04.33 


1643 


1.455 


0.519 


149.02 


-0.131 


0.037 


0.53 ± 0.03 


0.60 + 0.10 


0.43 + 0.04 


2.57 + 0.3 


4.2 ± 0.5 


84 


<4.2 


31365 


03-04.38 


03-04.40 


1653 


1.456 


0.520 


148.74 


-0.130 


0.036 


0.52 ± 0.03 


0.60 ± 0.10 


0.43 ± 0.04 


2.56 ± 0.3 


4.2 ± 0.5 


99 


<4.0 


31366 


03-04.44 


03-04.46 


1643 


1.456 


0.522 


148.46 


-0.130 


0.036 


0.52 ± 0.03 


0.59 ±0.10 


0.43 ± 0.04 


2.59 ± 0.3 


4.2 ± 0.5 


115 


<6.0 


31367 


03-04.51 


03-04.53 


1620 


1.456 


0.523 


148.18 


-0.128 


0.035 


0.52 ± 0.03 


0.59 ±0.10 


0.43 ± 0.04 


2.58 ± 0.3 


4.2 ± 0.5 


108 


<5.8 



° Two uvwl exposures were obtained for the observation 31332. 
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ground level. The peak of the radial profile is offset from the 
origin of the distance axis. This offset is described in more 
detail in Section 13.2.21 We extracted a background spectrum 
from the stacked event file, using a centrally positioned annu- 
lus with inner and outer radii of 500 and 645 arcseconds. We 
then extracted a spectrum using a centrally positioned circu- 
lar extraction region, with a radius of 349 arcseconds chosen 
to maximise the signal to noise. Both regions are shown over- 
laid on the exposure-corrected image of the stacked events of 
Figure [T] (upper panel). We applied the appropriate instrument 
XRT response matrix and effective area correction file (photon 
counting mode). 



0.5 



1 

Energy (keV) 



1.5 



Fig. 2. RASS backgrounds for the three sky positions of the 
comet (see text and Table 1), for 2009-01-28 (black), 2009-02- 
16 (red) and 2009-03-04 (green). 

and 31337 no exposures with the uvwl or v-band filters were 
obtained. 



3. Discussion 

3.1. Background and spectral analysis of the XRT data 

Comets are extended sources, often larger than the FOV of 
the observing instrumentation. Additionally, the X-ray back- 
ground can show considerable temporal and spatial variation. 
The background subtraction is therefore non-trivial, and we 
compare two different approaches here. First, we downloaded 
three rosat all-sky survey spectra (radius 1 degree) and the 
ROSAT PSPCc response matrix from the HEAS ARC background 
tooH where any known so-called Long Term Enhancements, 
due to time-v ariable SWCX (from e xospheric or heliosphere 
emission, e.g. Collier et al ]l2005l or iKoutroumpa et al.l l2006h 
have been removed dSnowden et al. Tl995h . The three point- 
ings were chosen to reflect the comet's position on the sky 
during the three observing periods (see Table 1). We plot the 
diree background spectra (from OBS JDs 4973969, 4917081, 
4814705) in Figure|2] The spectra show a considerable amount 
of variation. It was decided however, that creating three sepa- 
rate spectra for the comet would be infeasible due to the poor 
statistics that would be obtained. 

The second method, which we deemed more reliable, con- 
centrates instead on a spectrum obtained from the combined 
spectrum from all the available data for the comet. We show 
a radial profile of the exposure corrected image in Figure |3] 
Distances are quoted in arcseconds from the centre of the de- 
tector (which, following the procedure described in Section lZTI 
is the centre of the comet according to the ephemeris trans- 
formations that we had applied to the X-ray data). The radial 
profile allowed us to constrain the area in the detector plane 
that we considered to contain emission from the comet by se- 
lecting a radius at which the profile drops to the apparent back- 



0.0010 h 




200 300 400 

Distance (arcseconds) 



Fig. 3. Radial profile of the 0.3 to 1 keV combined XRT 
exposure-corrected image. The dashed vertical lines indicate 
the radius of the spectral (green) and the boundaries of the 
background extraction regions (red), chosen to maximise the 
signal to noise and to avoid subtracting signal coming from the 
comet. 



The SWCX emission process is characterised by line emis- 
sion in the soft X-ray band. Using the XSPECO spectral fitting 
package we fit an un-absorbed SWCX-like simple model to the 
background-corrected spectrum (shown in Figure 0, incorpo- 
rating three zero-width Gaussian lines, letting the energies and 
the line normalisations vary. The lines were found at 0.38 keV, 
0.52 ke V and 0.64 keV, which we attribute to a blend of helium- 
and hydrogen-like carbon and nitrogen lines, nitrogen and oxy- 
gen lines, and oxygen respectively. The reduced-Cash statis- 
tic (maximum likelihood-based statistic for Poisson data. ' CashI 
[1979) was 1.3 for 63 degrees of freedom with a not unreason- 
able goodness of fit of 35%. The goodness of fit should be ap- 
proximately 50% had the observed spectrum been produced by 
the model. The un-absorbed flux between 0.3 and IkeV was 
found to be 4.3 + 1.3 x lO"'-' ergscm"^ s"' (95% confidence). 
We also note the background-subtracted count rates in two en- 
ergy bands of interest: 4.1 + 1.1 xlO"^ ct s"'(0.3 to 0.5keV) 
and 4.1 + 1.6 xlO"^ ct s-i(0.5 to 0.7 keV). We look into the 
variation of flux with respect to time in Section l3.3.3l in com- 



^ http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg.pl 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html 
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parison to the behaviour of the solar wind at the time of the 3.2. Morphology. UV and X-ray 
observations. 

3.2.1 . Coma and tails 



0.08 
0.06 
0.04 
0.02 


0.02 



1 






1 ' 


1 ' 


D ^ 




H 












1 


1 



0.4 



0.6 0.8 
Energy (keV) 



Fig. 4. Upper panel: background subtracted spectrum from the 
combined comet Lulin XRT data set. The model fitted to the 
data, consisting of three zero-width Gaussian lines as described 
in the text, is shown by the solid line. Lower panel: residuals to 
the model fit. 



Several sophisticated methods have been developed to 
analyse cometary SWCX spectra (Beiersdorferetal. 2003[ 
Kharchenko et al. 2003 : Krasnopolskvl 20061 : Bodewits et al. 
2007h . We chose to apply a much simpler model here due to 



the low-quality spectrum achieved, even after stacking of the 
entire data set to improve the signal to noise. We do not observe 
any emission abo ve ~0.8keV, in contr ast to comets within the 
Chandra survey dBodewits et al. 2007 ). which consistently re- 
quired higher-order transitions of O viii plus two Ne lines at 
0.90 keV (Neix) and 1.024 keV (Nex). The Swift XRT has 
an effective area of ~65cm^ at ~0.8keV, much lower than 
the ~300cm^ at the same ener g y of t he Chandra ACIS S3 
chip used in the iBodewits et alj (120071) survey. However, the 
Swift XRT effective area increases to a maximum of ~ 120cm^ 
at about ~ 1.5 keV (compared to ~700cm^ for the Chandra 
ACIS S3 chip), implying that the lack of emission above 
this energy is due to incoming solar wind compositional ef- 
fects, which we discuss more in Section 13.3.31 Although the 
poor signal-to-noise ratio inhibits interpretation of the signal at 
higher energies, the lack of such a signal may suggest that the 
comet interacted with a solar wind with a charge state distri- 
bution reflecting a very low freeze-in temperature. The spectral 
shape of the model applied here mos t closely resembles that 
of comets 2P/ Encke (iLisse et al ]l2005h and 73P/Schwassmann- 
Wachmann 3 ( Wolk et al.l2009h of the Chandra sample, which 
sampled the coolest wind of that survey. We comment in more 
detail about the composition of the interacting solar wind in 
Section [333] 



UVOT images for every observing day are shown in Figure |5] 
The left most column shows the comet in the uvwl band for 
each observing periods (panels A, C and F). In this band the 
comet appears very symmetric around the nucleus, confirming 
that most of the emission seen in the uvwl band comes from 
OH molecules (Section r3.3.1l ). The comet's morphology in the 
v-band is driven by the changing viewing geometry. No v-band 
observations were obtained on 2009-01-28. On 2009-02-16 the 
dust coma is elongated in the solar direction. Subtracting the 
azimuthal profil e as determined by the meth od described in 
Section 122] (c.f. Schleicher & Famhamll2004 1 reveals both the 
dust tail (westward) and the anti-tail (eastward); the latter con- 
sists of long lived dust emitted near perihelion (panel E). Lulin 
moves in a retrograde orbit and Earth crossed the comet-Sun 
plane around 2009-02-26. On 2009-03-04, the observing ge- 
ometry therefore no longer allows easy separation of particles 
along different trajectories, resulting in a dust tail in the anti- 
solar direction (shown by the elongated contours in panel G). 

On 2009-01-28, a well defined, narrow tail can be seen to 
extend westwards (best seen in panel B, where the image is 
divided by the azimuthal profile). It is not visible in the 2009- 
02-16 images (panel E). The nature of this narrow tail was first 
discussed in Bodewits et al. (201 1). In this paper it was argued 
that from the geometry the narrow tail could be either the dust 
tail or the ion tail (both clearly visible in wide field amateur 
images obtained in this period), but that while the uvwl filter 
has a significant red-leak, it is not very sensitive to the emission 
usually associated with the blue CO^ comet tail system (4000 
- 5000 A). It was concluded that it was more likely that the 
narrow feature arose from continuum emission from dust or 
ice. No asymmetric feature, however, is visible in the uvwl and 
v-band azimuthal profile-subtracted images on 2009-02-16. 

As we argued above, the UV filter transmission means 
that the narrow tail is likely to be small-grained dust rather 
than ions. The geometric conditions and the fact that the tail 
was detected in the UV observations in January but not in 
February suggests that it is the result of enhanced activity of 
these grains around perihelion. The dust, released around per- 
ihelion and pushed away from the Sun, was projected to the 
west in January, matching the direction of the observed tail. By 
the time of the February observations, the Earth had moved to 
the point that we were looking almost straight down this tail. 
Unfortunately, the lack of UVOT v-band observations preclude 
the comparison of larger dust grains between the two dates. 

3.2.2. Solar wind interaction in X-rays 

As discussed in the previous section, the OH distribu- 
tion in the coma was very symmetric and homogeneous 
during our observations. Past modelling has indeed indi- 
cated a transition between thick and thin regimes between 
10^^ and lO^^mol. s"'. Lulin's large gas production rate (6- 
8 xlO^**mol. s"^ Sec. 13.3.1b made it collisionally thick to 
charge exchange by solar wind ions (c.f. iLisse et al.ll2005l and 




uvwl 
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Fig. 5. UVOT observations of Comet Lulin. A. 2009-01-28,uvwl. B: image A divided by coma profile. C: 2009-02-16, uvwl. D: 
2009-02-16, v-band. E: image D with the coma profile subtracted. F: 2009-03-04, uvwl G: 2009-03-04, v-band. All images are 
oriented in the same way (east left, north up) and have the same angular scale (UVOT's FOV is 17 x 17 arcminute). Iso-intensity 
contours are shown on a logarithmic scale to enhance fainter features. In panel F the solar direction vector and the orbital vector 
are coincident. 



Bodewits et al.ll2007h . It is therefore to be expected that when 
viewed via X-ray emission, the comet would appear as a cres- 
cent when observed at phase angles of around 90 degrees, 
with a significant sunward offset between the maximum of the 
observed X-ray emission and the nucleus of the comet (cf . 
Hyakutake (Lisse et al. 1 993); Linear S4 dLisse et all 1200 ll): 
C/ 2001 WMl dWegmann & Dennerlll2005h V From the models 
by Bodewits et aP 20071 most charge exchange reactions in the 



coma of a comet with a gas production rate of ~ 10^^ mol. s"' 
would occur lO^* - 10^ km from the nucleus (at a distance of 
1 AU from the Sun). 

In Figure [1] (lower panel) we plot X-ray contours over a 
UVOT UV image of the comet (from the first observation, 
31332 using filter uvwl). As this X-ray image is constructed 
from images taken at different Sun-Earth-comet angles and 
geocentric distances (resulting in diff'ering angular extents of 
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the expected X-ray emission), caution should be exercised 
when interpreting the implied morphology. A very rough cres- 
cent shape can be seen in Figure [T] (lower panel). Assuming 
the exposure-weighted distance of 0.685 AU to the comet we 
find that the X-ray maximum is offset in the sunward direc- 
tion from the UV maximum by ~35,000kin (foreshortened 
by the observing geometry ). While in line with other co met 
observations (iWegmann & Dennerj2005 : Dennerl et al .l2003h . 
this stand off" distance is greater than expected by the model 
results of iBodewits et al.l (12007.) . who predict the brightness 
peak around 10,000 km for Q = 10^^ molecules/s at a helio- 
centric distance 1 AU, and less than predicted by the model of 
(WegmanrLet aL. (2004.) . who predict 66,000 km under the same 
circumstances. 



3.3. Temporal variations 

Variations in the comet's X-ray luminosity are driven by vari- 
ations in the comet's gas production rate, in the solar wind ion 
flux, and by the solar wind ion content. Along with its X-ray 
observations. Swift simultaneously obtained UV/optical obser- 
vations of the comet, which provide a measure of the comet's 
gaseous activity. The solar wind is currently sampled by vari- 
ous spacecraft and data from these spacecraft are the best prox- 
ies for near-Earth comet environments. In this section we in- 
terpret our observations in terms of the variability of both the 
cometary neutral gas content and the solar wind ion flux. 

3.3.1 . Gas production rates 

At the wavelengths covered by UVOT, comets are seen in sun- 
light reflected by cometary dust, with several bright molecu- 
lar emission bands superposed. The uvwl filter is well placed 
to observe the three very strong OH vibrational transitions be- 
tween 280 and 312nm. Lulin' s spectrum was only sl i ghtly red- 
dened (< 10% per 1000 A: iBodewits et al.1 1201 ll iLin et al 
2007b . We therefore assumed an un-reddened solar continuum, 
and convolved this with the uvwl and v-band filter transmis- 
sions to determine how much the continuum contributes to the 
uvwl flux. If we further assume that the flux in the v-band filter 
is dominated by continuum emission (while in truth it is con- 
taminated by the fluorescent emission of various molecules, 
predominantly C2 and NH2). The OH flux is then given by 
Equation [H 



^'0ff = ^'m.vvl -(0.15X^1/) 



(1) 



where Fu^ki and Fy are the fluxes measured in the uvwl and v- 
band filters. The factor 0.15 originates from the ratio Fumi/Fy 
for a solar spectrum. From Table 1, on 2009-02-16 the contin- 
uum contribution (here assumed to be any contributions other 
than from OH) to the flux measured in the uvwl band was 
approximately 21%; on 2009-03-04 it was 17%. We consider 
the possibility that there is a non-negligible contribution from 
the fluorescence emission of C2 to the v-band measured flux. 
Using the total C2 flux wit hin this band determi ned from our 
work with the Swift grism (IBodewits et al. 201 lb . we estimate 
the contribution to be 0.25 W m~^ sr~' . This corresponds to an 



over-subtraction of 8% in the calculation of the OH flux using 
Equation[T] As this is well within the errors quoted in Table 1, 
we consider this contamination of marginal consequence. 
Using the molecules' fluorescent rate 



dSchleicher & A'Heamlll988L the number of OH molecules in 
the aperture can be derived (IFeldman et al. I l2005h . The results 
are listed in Table 1 . To derive production rates, we modelled 
the OH and water distribution, and compared this with the 
measured OH number densities. This model is based on the 
Haser model (Haser 1957; Festou 1981) and is discussed in 
greater detail in IBodewits et al.l (|201 lb . 

The use of broad band filters implies that our results 
are a relatively crude estimate of the comet's water pro- 
duction and we estimate the systematic error to be around 
25%. The resulting OH production rates are summarised in 
Table 1. No v-band observations where obtained on 2009- 
01-28, and without accounting for the continuum removal 
we find a 3-sigma upper limit of Qoh < 1 xlO^^mol. s"\ in 
good agreement with the production rates of 5.8+0.6 (1- 
sigma) to 6.9+0.6x10 ^^ mol. s ' measured with the grism 



d Bodewits etall l201lb . On 2009-02-16 and 



2009-03-04 
and 



the OH production rates were 5.5 
4.2+1.0x10^^ mol. s"', respectively. 

The OH production rates were converted to water produc- 
tion rates assuming a branching ratio for t he formation of OH 



from H2O of 85% dHuebner et al.ll 19921; ICombi et al.ll2004b . 



In Figure |6] we show these rates versus time, for this work 
and other measurem e nts (C ombi et al. 2009; Bone v et al .I2OO9I 



Ootsubo etaPbOlOl: IBodewits et al.ll201lL and D. Schleicher 



(private communication). Our measurements indicate that the 
comet's activity was slowly decreasing after it had reached per- 
ihelion on UT Jan. 10.64, 2009, in good agreement with other 
measurements. 

Several studies found that gas production rates varied on a 
day-to-day basis by a s much as 50%, suggesting strong diurnal 
eifects on the comet ([Bodewits et alj|201 It ICombi et al. I I2OO9I; 



Knight & Schleichej l2009h . Based on the phasing of CN jets. 



Knight & Schleichei (l2009b found a rotation period for comet 



Lulin of 42+0.5 hr Even though our UVOT broadband obser- 
vations on 2009-02-16 and 2009-03-04 sampled the comet ev- 
ery hour for about half a day and thus we did not observe the 
comet for a full rotation period, our data showed no significant 
short term flux variations. 

To estimate the comet's neutral gas contribution to its X- 
ray variability one needs to know the number of neutral coma 
molecules in the FOV It can be assumed that the X-rays are 



mostly driven by SWCX with H2O and OH ( IBodewits et al 
12OO61) . We therefore need to know the sum of the number of 
H2O and OH molecules in the aperture. If the aperture were 
large enough, the total number of OH molecules would be 0.85 
times the number of H2O molecules (the branching ratio of the 
H2O to OH photodissociation process), and the total number of 
neutral molecules in the FOV would thus be 2.2 times the num- 
ber of OH molecules present. However, in apertures smaller 
or equivalent to the Haser scale-lengths for H2O to OH dis- 
sociation (of the order of 10^ km at 1 AU from the Sun), the 
number of H2O molecules will be larger than the number of 
OH molecules. For the apertures used here (~ 3.1 x 10^ km. 
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Fig. 6. Gas production rates of H2O vs time, calculated for 
our observing periods on 2009-02-16 and 2009-03-04 (solid, 
black diam onds). We also plot H2O production rates from 
Combi et alJl2009l (circl es), Schleicher 2009 ( private commu- 
nication, open squar es), Bodewits et al. 20 ll (filled squares) 



Bonev et al. 2009 



Ootsuboet^ I2OIOI (filled triangle) and 
(solid, grey diamond). The dashed vertical line indicates the 
date of the perihelion passage of the comet. 

~ 1.45 X 10^ km and ~ 1.51 x 10^ km radius, for the three ob- 
serving periods), the total number of H2O plus OH molecules, 
calculated by Haser modelling, are 1.54, 1.56, and 1.57 times 
the number of observed OH molecules on 2009-01-28, 2009- 
02-16, and 2009-03-04, respectively. This factor incorporates 
the lifetimes of the H2O and OH molecules with cometary dis- 
tance from the Sun and accounts for the proportion of H2O 
compared to OH observed depending on aperture size. For the 
2009-01-28 measurement the aperture is much larger than the 
Haser scale length, whereas the factors for 2009-02-16 and 
2009-03-04 are adjusted to account for the smaller apertures 
used. 

To estimate the number of water and OH molecules on 
2009-01-28, when we did not have any v-band images to sub- 
tract the continuum emission from the uvwl filter, we used the 
OH production rates derived from the grism observations and 
used our Haser model to estimate the number of molecules in a 
400 arcseconds radius aperture. 

3.3.2. Solar wind 

We took data from stereo a to investigate the behaviour of the 
solar wind in early 2009, shown in Figure |7] We plot the so- 
lar wind proton speed and flux, mapped to the approximate lo- 
cation of the comet, stereo a is found ahead of the Earth in 
its orbit on the ecliptic plane. In order to map the solar wind 
from STEREO A to the position of the comet at the three ob- 
ser ving period s , we used the time shift pro c edure described 
by iLisse et alj d 19991) and iNeugebauer et all (|2000|) . The cal- 




40 50 
Day of year 200S 




40 50 
Day of year 2009 

Fig. 7. Solar wind parameters, taken from stereo a, that have 
been co-rotationally mapped to the location of the comet. 
Upper panel: proton velocity. Bottom panel: proton flux. The 
three observing periods are highlighted within the light-grey 
shaded areas. 



culations are based on the comet ephemeris, the location of the 
comet and stereo a, and the bulk solar wind velocity measured 
by STEREO A. With this procedure, the time delay between an el- 
ement of the co-rotating solar wind arriving at the proton mon- 
itor and the comet can be predicted. Although this mapping is 
an approximation that does not take propagating shocks or lat- 
itudinal structures in the wi nd into account (Neugebaue r et al.l 
2OOOI: iBodewits et al.ll2007h . the relatively close proximity of 
the comet and its low heliocentric latitude (<1 degree during 
our observations) aid its reliability. 

On 2009-01-28, there is a delay of 2.5 days between the so- 
lar wind measured by stereo a and the solar wind arrival at the 
comet, mostly due to the longitudinal difference. On 2009-02- 
16, this delay is 0.8 days, and on 2009-03-04 the comet leads 
STEREO a by 1.8 days. It is of note that the uncertainty in the so- 
lar wind arrival time is the largest in January, when the comet 
was still at 1 AU from Earth versus 0.5 AU in March. 

The averaged solar wind proton fluxes, were 15, 5 and 14 
X 1 0^ cm-2 s- 1 for 2009-0 1-28, 2009-02- 1 6, and 2009-03-04 re- 



spectively. The average solar wind proton speeds for the three 
observing periods were ~ 370, 280 and 420km s"', implying 
solar wind densities of 4.2, 1.9 and 3.4 cm"^. Inspection of long 
term solar wind data indicate the presence of several CIRs (a 
compression region followed by a rarefaction, as described in 
Section [1]) during the early part of 2009. The most notably of 
these occurred just after the second observing period, as indi- 
cated by sharp increases in solar wind proton velocity. In addi- 
tion a weaker CIR occurred during the third observing period, 
shown in the upper panel of Figure |7] The low solar wind pro- 
ton densities may be explained by the presence of these CIRs. 
The solar wind data further show that a CIR reached the comet 
around 2009-02-05 (day of year 36). Considering the uncer- 
tainty in the time delay of our propagation model, we consider 
it likely that it was this CIR that caused a disconnection event 
at the comet, reported bv lShi et al ] |2011[ 
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3.3.3. X-ray lightcurve 

In Figure [8] (upper panel) we plot the background-corrected 
XRT fluxes in the energy band 0.3 to 1 keV for each observa- 
tion (accounting for the areas of the spectral extraction regions 
to convert to the stated flux units), or the upper limit where 
appropriate. We also show the predicted X-ray flux, calculated 
from the product of solar wind proton flux and the number of 
neutral molecules in the FOV, corrected for the Earth-comet 
distance (Equation|2]i: 



A^ga.s * f sw(rh) 

A2 



(2) 



A^aas is the total number 



where A; is a normalisation constant, 
of water group molecules available to contribute to the charge 
exchange process derived from our UVOT observations (see 
Section [3.3.1l ). Fsw('"h) are the solar wind proton fluxes at the 
comet (as described in Section I3.3.2l i. and A are the Earth- 
comet distances (as listed in Table 1). We normalised to the 
flux measured in the tenth XRT observation, as this data point 
had both a constrained X-ray and UVOT measured flux. In 
Figure [8] (middle panel) we plot the solar wind proton flux, 
at the three observing periods as described below. The lower 
panel of Figure |8] shows A^gas over the periods of the observa- 
tions. 

There are several competing processes that contribute to the 
predicted X-ray flux as given in Equation|2] The gas production 
rate was 50% larger in January compared to the other observing 
periods, however the comet was at half the distance from Earth 
in February and March. The average proton flux was compara- 
ble during 2009-01-28 and 2009-03-04, and about 40% lower 
during our 2009-02-16 observations. The contributing parame- 
ters to the predicted lightcurve compete to cancel each other 
out, resulting in a relatively flat predicted X-ray lightcurve. 
Thus our measurements do not exhibit any large variations, and 
our predicted and measured lightcurves are in reasonably good 
agreement. 

The solar proton flux is only a proxy for the heavy ion flux 
of the solar wind and this heavy ion content may vary con- 
siderably due to variations in either solar wind atomic abun- 
dance or ion temperature and charge state. Thus we looked 
at additional available solar wind O^^ flux data during the 
time of our obs ervations (e.g. from the Advanced Composition 
Explored (ace), IStone et alJI 19981) . located at Lagrangian Point 
LI). Unfortunately the available data were sparse and of low 
quaUty. 

Had the comet interacted with a ICME (a large cloud of 
plasma ejected from the Sun, as described in Section [TJ during 
our observations, we would have in general expected a large 
increase in the measured X-ray flux and spectral hardness due 
to the presence of highly charged ions (as shown by the exam- 
ples of an ICME interacting with th e exo sphere of the Earth; 



Ezoe et alj boi it ICarter et al.1 12OIOI) . Approximately 10% of 



ICM Es appear to exhibit only weak compositional anoma- 
lies dRichardson & Cane 2004). The long-term trend data of 
Figure [T] suggested that during at least one, if not two of the 
observing periods, the comet had sampled a CIR. CIRs occur 
when a fast solar wind stream piles up against a slow solar wind 
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Fig. 8. Upper panel: XRT 0.3 to 1 keV background-corrected 
fluxes (diamonds) or upper limits (downward arrows) where 
appropriate, plotted against a time-axis given in the day of year. 
We also plot the predicted X-ray flux (red crosses, normalised 
to the tenth observation) as described in the text. Middle panel: 
solar wind proton flux as measured by stereo a for early 2009. 
Lower panel: Ngas (open diamonds), or upper limits (downward 
arrows) where appropriate. 



stream. CIRs show average compo sitional sig natures similar to 
average fast and cool solar wind dMason et al. 2008). Hence, 
the comet may have sampled, or partially sampled, a fast, cool 
plasma that would have exhibited lower abundances of highly- 
charged heavy ions that are responsible for X-ray emitting, 
charge-exchange interactions. The X-ray spectrum of Lulin is 
indicative of a higher solar wind charge state than the very sof t 



X-ray spectrum of comet 17P/Holmes ( Christian et al. 20101) 



which showed no X-ray emission above 0.4 keV and was pos- 
tulated to be the result of an interaction with the cool, polar 
solar wind. 
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4. Conclusions 

We have used the unique capabilities of the Swift observatory to 
simultaneously observe comet C/2007 N3 (Luhn) via UV and 
X-ray emission. We have carefully transformed X-ray events 
from 14 Swift XRT pointings towards comet Luhn to produce 
the optimum exposure-corrected X-ray image. We have also 
analysed Swift UVOT exposures when available and used this 
data to observe variations in the morphology of the comet in 
particular with regards to the dust and ion tails. We see the peak 
of the cometary X-ray emission offset from the peak in the UV 
emission, and this offset is displaced towards the Sun, indicat- 
ing that the coma was coUisionally thick to charge exchange 
from the nucleus out to about 35,000km, as expected for a 
comet with a gas production rate of QhjO = 6 xlO^^ mol. s"' . 

A model applied to a background-corrected X-ray spectrum 
was consistent with emission from charge exchange, although 
individual lines were heavily blended. The spectrum, only dis- 
cernible below 1 keV, also suggested that the comet had sam- 
pled a cool solar wind, with less highly ionised oxygen than a 
quiescent (warm) equatorial solar wind. Such a spectrum may 
be expected if the comet had encountered a stream of fast wind 
originating at the cold base of the corona and flowing through 
a coronal hole. We have used the solar wind proton flux as a 
proxy for solar wind activity throughout this work, as investi- 
gations regarding solar wind minor ion fluxes were hampered 
by a lack of high cadence data from spacecraft solar wind ion 
monitors during the Swift observations. What heavy ion data 
trends that do exist suggest there were indeed numerous CIRs 
during the period of the Swift Lulin observations. 

We combined our simultaneous X-ray and UV observa- 
tions with measured solar wind data to explain the measured 
X-ray lightcurve. We use data from the Swift UVOT, employ- 
ing the uvwl and v-band filters to measure the number of 
OH molecules in the FOV and estimate the number of water 
molecules released by the comet. We used the water production 
rates of the comet and solar wind data mapped to the location of 
the comet to estimate the expected X-ray flux. We demonstrate 
that Lulin's X-ray brightness was determined by the comet's 
activity and the observing geometry, as weU as variations in 
the solar wind. 
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